An Improved Prediction Model of Vortex Shedding Noise from Blades of Fans 
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The main source of the noise of an axial flow fan is the fluctuating pressure field on blade surfaces caused by the 
shedding of vortices at the trailing edge of blades. An analytical model to predict the vortex shedding noise gen- 
erated at the trailing edge of blades of axial flow fans was proposed by Lee in 1993. In this model, for mathemat- 
ical convenience, an idealized vortex street is considered. However, the agreement between the analytical results 
and the experimental data needs to be improved because of the simplification about the Karman vortex street in 
the wake of blade. In the present study, a modified model is proposed based on the prediction model by Lee. The 
boundary layer theory is used to analyze and calculate the boundary layer development on both the pressure and 
the suction sides of blades. Considering the effect of boundary layer separation on the location of noise source, 
the predicted overall sound pressure level compares favorably with the experimental data of an axial fan. In the 
calculation of A-weighted sound pressure level (La), considering the effect of static pressure on radiate energy, 


the predicted broadband noise with the modified model compares favorably with the experimental data of a mul- 


tiblade centrifugal fan. 
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Introduction 


Considering the environments of working and life, the 
problem of noise generated by fans has been obtained 
great concern over the recent years. The aerodynamic 
noise of fans is very dominant in fan noise, rather than 
the mechanical noise caused by vibration and the elec- 
tromagnetic noise caused by motor [1]. Lighthill [2] 
showed that the vortex shedding noise generated from the 
wake of blades of fans is the main source of noise. 

Several prediction models of noise based on Karman 
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vortex street have been widely utilized for the prediction 
of noise from axial fans or centrifugal fans. Sharland [3] 
proposed that the broadband noise of an axial flow fan 
can be described by three mechanisms, surface pressure 
field arising from the turbulent boundary layer, influence 
of the vortex shedding from the surface of a rigid body, 
and the incoming flow with initially turbulent flow. Lee 
[4] proposed an analytical model for predicting the vor- 
tex shedding noise generated from the wake of axial flow 
fan blades. Wang [5] developed a noise model and opti- 
mized the algorithm of boundary layer development. 
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Wang [6] presented an analytical model for predicting the 
noise of plenum fan based on the structure characteristics 
of centrifugal impeller. Khelladi [7] showed that in a far 
field a monopolar source is equivalent to a dipolar source 
induced by a uniform distribution of the load on the en- 
tire moving surface. Sasaki [8] proposed a prediction 
theory for the broadband noise generated from a multib- 
lade fan based on the characteristics of the Karman vor- 
tex street. Capece [9] investigated the unsteady flow 
phenomena resulting in the differences in the instantane- 
ous rotor blade wake data. In these researches, it is well 
agreed that the noise generated by vortex shedding is the 
primary source of noise for small fans. 

In the present work, an improved predicting equation 
of broadband noise generated by fans is proposed. With 
this improved equation, the predicted broadband noise 
compares favorably with the experimental data of a mul- 
tiblade fan in the literature [8]. 


Model and Method 


Lee developed an analytical model for predicting the 
vortex shedding noise generated from the wake of axial 
flow fan blades [4]. 

According to the stability analysis of the Karman vor- 
tex street, the distance between vortex rows can be esti- 
mated by [4], 

b=0.66,, + 0.66,, +d (1) 


where b is the distance between vortex rows, d is the 
trailing edge thickness, ôs and p are the boundary layer 
thickness on suction surface and pressure surface, re- 
spectively. 

s= 2 =0.281 (2) 
a 

where s is a constant, and a is the distance between 
vortices. 

The vortex strength can be written as [4], 


V ‘sa - (say —2na(0,, + Op (27s tanh(zs) — D| 


27s tanh(zs) —-1 
(3) 

where K is the strength of each individual vortex, V 
represents the main flow velocity, and 6,, and 6,, are the 
momentum thickness of boundary layers on suction sur- 
face and pressure surface, respectively [4]. 

The whole vortex street has an advance velocity, V; = 
(K/2a) tanh (as), then the vortex shedding frequency is 
given by [4], 

f=V,-V)la (4) 

The Strouhal number based on the trailing edge thick- 
ness is written as [4], 

St = fd/V, (5) 


The pressure distribution can be written as [4], 
P(x,t) = 
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where p is the density of air, c is the chord length, and 
ois the angular velocity, 
o=2n(V,-V;)/a (7) 
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In Eq. (6), &, is the distance from the centre of vortices 
to the origin which is located at the centre of blade. 

In Lee’s study, for mathematical convenience, an idea- 
lized vortex street (see Fig. 1) has been considered, 
where c is the chord length. In Fig. 1, the origin o is set at 
the centre of airfoil. The noise source is considered to be 
at the trailing edge of blade since it is assumed that the 
vortex shedding occurs at the trailing edge. 

However, since the boundary layer flow at the suction 
side of the blade is under adverse pressure gradient at an 
attack angle, it is easy to be separated from the blade 
surface (see Fig. 2). When the boundary layer is sepa- 
rated, vortices would be generated immediately from the 
separation point. At high Reynolds number, vortex shed- 
ding could be produced upstream at the separation point 
rather than at the trailing edge of blades. 

Figure 3 shows the flow around an airfoil which is 
taken from a cascade and is studied at the attack angle 8 
of 17.5 degree [10]. 


Fig. 1 Idealized model of vortex street 


Fig. 2. Modified model of vortex street 


Fig. 3 Schematic diagram for airfoil 


As discussed above, considering the separation of 
boundary layer, the location of noise source may be 
moved from trailing edge to the separation point. Taking 
into account of this phenomenon, an improved calcula- 
tion equation of noise is proposed as follow, 

É, =Xstna/2 (9) 

In Eq. (9), the first term x, stands for the distance from 
the origin o (the center of blade) to the separation point S, 
and the second term represents the distance between a 
shedding vortex and the separation point S. 

The sound power e radiated from the incoming flow 
field can be expressed as [4], 
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where agis the speed of sound, the constants by and b, 
are 0.19 and 0.71, respectively [11], and the first 50 terms 
taken from the series can meet the actual computing re- 
quirements [6]. 

In reference [8], Sasaki et al. did experiment for a 
multiblade centrifugal fan and compared the experimen- 
tal result with the prediction with Lee’s classic model [4]. 
It is found that the agreement between the theory and the 
experiment is very poor. Therefore, an improvement is 


(10) 


needed on Lee’s theory for prediction of broadband noise. 


In the following, considering the effect of the static pres- 
sure in vortex flow domain, the equation for calculating 
the frequency of broadband noise generated by vortex 
shedding is proposed. The frequency of vortex shedding 
can be written as [8], 
f =StxV,d" (11) 
where St expresses the Strouhal number and the V, is 
the relative velocity at impeller outlet. 

The total pressure represents the total mechanical en- 
ergy of fluid. The level of total pressure at impeller outlet 
can be used to characterize the noise level. The total 
pressure consists of the static pressure and the dynamic 
pressure, 
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where P, is the total pressure, P,, is the static pressure, 
and V, is the absolute velocity at impeller outlet. 

In the flow model [8], the flow has been divided into 
two domains. One is a domain with the vortex flow at the 
bellmouth side (vortex flow domain), and the other is the 
mainstream zone where the outflow is close to the hub 
side (mainstream domain) [8]. The static pressure in vor- 
tex flow domain is clearly larger than that in main flow 
domain from the experiment in [8]. 

In Lee’s model [4], only the dynamic pressure is in- 
volved in the calculation of radiate energy, and the effect 
of static pressure is neglected. In order to considerate the 
influence of static pressure, a pseudo-velocity Vis pro- 
posed in this study by dimensional analytical method, 
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which is used to correct the main flow velocity. 

It is assumed that the influence of static pressure on 
the overall sound pressure level only exists in vortex 
flow domain. Then the critical value of the static pressure 
used to calculate the pseudo-velocity V” is the minimum 
value in vortex flow domain. On the other hand, it is as- 
sumed that the influence of static pressure to the overall 
sound pressure level is negligible in mainstream domain. 
Thus, the critical value about static pressure is also the 
maximum value of static pressure in mainstream domain, 


p(en 
p 
where P, is the critical pressure, which is taken as the 
minimum value of static pressure in vortex flow domain. 
As such, the corrected velocity is defined by, 
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Then the total radiated energy can be expressed as, 
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There is a relationship between the sound power and 
the sound pressure p [4], 
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where / is the distance from the test point to the noise 
source. 
The Eq. (18) can be re-written as [6], 
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where B is the number of blades. 
The overall sound pressure level is defined as, 


2 
SPL = 201g = 101g", (20) 
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where p,ep is a reference pressure which is equal to 
2x10” Pa. 


Results and Discussions 


Fukano et al. did experiment for an airfoil-shaped axi- 
al fan [10], and compared the predicted overall sound 
pressure with Lee’s model and the experimental data. The 
result shows that the predicted overall sound pressure is 
lower than that from experiment [10], as shown in Fig. 4. 

The reason for the discrepancy between the theory and 
the experiment may be due to the separation of boundary 
layer. The prediction result with the improved model 
which considering the effect of boundary layer separation 
on the suction surface is also plotted in Fig. 4. It is seen 
that the modified model agrees better than Lee’s original 
model with the experiment. Thus, considering the effect 
of boundary layer separation, the sound pressure level 
can be predicted more accurately. 
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Fig. 4 Overall sound pressure of the axial fan versus number 
of blades; The experimental data is taken from [10] 


Sasaki et al. did experiment for a multiblade centri- 
fugal fan [8]. In their study, the flow field in the fan is 
simulated with commercial software Fluent. The simula- 
tion results of the efficiency and total pressure of the fan 
are compared with the experimental data as shown in Fig. 
5 and Fig. 6, respectively. In Fig. 5, the ordinate is the 


efficiency y and the abscissa is the flow rate coefficient ø. 


In Fig. 6, the ordinate is the total pressure P, and the ab- 
scissa is the flow rate coefficient ø. It can be seen that the 
simulated results agree well with the experimental data. 
Figure 7 shows the distribution of relative velocity and 
absolute velocity versus the impeller width at the trailing 
edge of blade (where z=0 corresponds to the bellmouth 


side and z=1 corresponds to the hub side). As mentioned, 
the flow has been divided into two domains, main do- 
main (hub side) and vortex domain (shroud side). It can 
be seen from Fig. 7 that the velocities in main flow zone 
is higher than those in vortex zone. This trend is in 
agreement with the experimental data [8]. 
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Fig. 6 Measured and simulated total pressure versus flow rate 
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Fig. 7 Distribution of velocity along the blade height (V, and 
V, expresses the absolute and relative velocities, re- 
spectively) 


If the vortex shedding noise generated from the wake 
of fan blades is the main noise source, the vortex shed- 
ding frequency is also the frequency of broadband noise. 
Thus, the relative velocity is used to calculate the fre- 
quency of broadband noise in Eq. (11). In the following 
calculation, the absolute velocity will be used to calculate 
the radiate energy in Eq. (17) because of the radiate 
energy is closely related to the absolute velocity. 

Figure 8 is the distribution of static pressure at the 
impeller outlet along the trailing edge of blade (where 
z=0 corresponds to the bellmouth side and z=1 corres- 
ponds to the hub side). It can be seen that the static pres- 
sure in vortex flow domain is clearly higher than that in 
main flow domain. In the actual calculation of radiate 
energy, only the effect of static pressure in vortex zone is 
considered, while the effect of static pressure in main 
stream zone is not taken into account. 
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Fig. 8 Distribution of static pressure along the blade height 


In Fig. 8, it can be seen that the static pressure is rela- 
tively high when 0<z<0.4. Taking into account the veloc- 
ity along the trailing edge of blade, the region of 0<z<0.4 
is regarded as the vortex flow, and the region of 0.4<z<1 
is thought to be the main flow region. This partition is 
consistent with the result of the literature [8]. Thus, the 
static pressure at z=0.4 is taken as the critical pressure, P,, 
in Eq. (15) for the detailed calculation of radiate energy. 

In Fig. 9, the predicted broadband noise with the im- 
proved model is compared with experimental data in [8]. 
Result predicted with Lee’s original model [4] is also 
included. The experiment on the fan was performed un- 
der the condition of 1400 r/min [8]. Thus, the discrete 
noise is generated in the blade path frequency about 
2333Hz. In addition, because of the noise in high fre- 
quency is generated by harmonic mainly, this model 
mainly pays attention to the broadband noise which is 
generated from 100 to 2000Hz. 

It can be seen from Fig. 9 that the frequency of 
broadband noise calculated by Lee’s model is deviated 
from the broadband noise area. 
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Fig. 9 Comparison of the predicted broadband noise with 
experimental data in [8] 


With the improved model, the predicted A-weighted 
sound pressure level (La) is plotted in Fig. 9, where the 
frequency is calculated with the Strouhal number by Eq. 
(11). In the present calculation, St=0.2 is employed ac- 
cording to reference [8]. Without considering effect of 
static pressure, the predicted result is expressed as 
“present V”. With considering effect of static pressure, 
the predicted result is expressed as “present V*” in Fig. 9. 

It can be seen from Fig. 9 that both calculation results 
in this study agree well with the experimental data in 
range of 100Hz~2000HZ. With considering effect of 
static pressure, the predicted result of A-weighted sound 
pressure level (La) favorably accords with the experiment. 


Conclusions 


In the prediction of noise of fans, Lee’s classic model 
is widely used by several researchers. However, this 
model compares poorly with experiments in some cases. 
In the present study, improvement on this model is car- 
ried out by considering effect of some factors. 

Firstly, the boundary layer theory is used to analyze 
and calculate the boundary layer development on both 
the pressure and suction sides of blades. Considering the 
effect of boundary layer separation on the location of 
noise source, a simple improved model is constructed. 
The predicted overall sound pressure level compares fa- 
vorably with the experimental data of an axial fan. 

Secondly, in the calculation of A-weighted sound 
pressure level (La), considering the effect of static pres- 
sure on radiate energy, an improved model is proposed. 
The predicted broadband noise with the modified model 
compares favorably with the experimental data of a mul- 
tiblade centrifugal fan. 


Acknowledgement 


This work is supported by National Natural Science 


Foundation of China (51206149, 51376055) and Zhe- 
jiang Province Science and Technology Innovation Team 
Project (2013TD18). 


References 


[1] 


[2] 


[5] 


Chen, X. D.: Centrifugal Fan, Chemical Industry Press, 
Beijing, pp.117—119, (2006). (in Chinese) 

Lighthill, M. J.: On Sound Generated Aerodynamically. I. 
General Theory, Proceedings of the Royal Society A, 
vol.211, pp.564—587, (1952). 

Sharland, I. J.: Sources of Noise in Axial Flow Fans, 
Journal of Sound and Vibration, vol.l, pp, 302-322. 
(1964). 

Lee, C., Chung, M. K., Kim, Y. H.: A Prediction Model 
for the Vortex Shedding Noise from the Wake of an Air- 
foil or Axial Flow Fan Blades, Journal of Sound and Vi- 
bration, vol.164, pp.327—336, (1993). 

Wang, D. B.: A Numerical Prediction Model for the 
Noise Generated from Low-pressure Axial Fans, Fluid 
Machinery, vol.29, pp.25—27, (2001). (in Chinese) 


[6] 


[7] 


[8] 


(11) 


Wang, W. X., Li, S.: Huang, D. T., Zhu, Z. C., Zhu, K. Q.: 
Numerical Prediction of Vortex Shedding Noise Radiated 
by a Centrifugal Impeller, Engineering Mechanics, vol.25, 
pp.217-222, (2006). (in Chinese) 

Khelladi, S., Kouidri, S., Rey, R.: Isom's Thickness Noise 
for Axial and Centrifugal Subsonic Fans, Journal of 
Sound and Vibration, vol.313, pp.1—6, (2008). 

Sasaki, S., Kodama, Y., Hayashi, H.: Influence of the 
Karman Vortex Street on the Broadband Noise Generated 
from a Multiblade Fan, Journal of Thermal Science, 
vol.14, pp.198-205, (2005). 

Capece, V. R., Fleeter, S.: Measurement and analysis of 
unsteady flow structures in rotor blade wakes, Experi- 
ments in Fluids, vol.7, pp.61—67, (1989). 

Fukano, T., Kodama, Y., Takamatsu, Y.: Noise Generated 
by Low Pressure Axial Flow Fans. II. Effects of Number 
of Blades, Chord Length and Camber of Blade, Journal of 
Sound and Vibration, vol.50, pp.75—88, (1977). 

Brooks, T. F., Hodgson, T. H.: Trailing Edge Noise Pre- 
diction from Measured Surface Pressures, Journal of 
Sound and Vibration, vol.78, pp.69-117, (1981). 


